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Abstract 

Ticks are one of the most important blood-sucking vectors for infectious microorganisms in humans and animals. When 
feeding they inject saliva, containing microbes, into the host to facilitate the uptake of blood. An understanding of the 
microbial populations within their salivary glands would provide a valuable insight when evaluating the vectorial capacity of 
ticks. Three tick species (Ixodes ovatus, I. persulcatus and Haemaphysalis flava) were collected in Shizuoka Prefecture of Japan 
between 2008 and 2011. Each tick was dissected and the salivary glands removed. Bacterial communities in each salivary 
gland were characterized by 16S amplicon pyrosequencing using a 454 GS-Junior Next Generation Sequencer. The 
Ribosomal Database Project (RDP) Classifier was used to classify sequence reads at the genus level. The composition of the 
microbial populations of each tick species were assessed by principal component analysis (PCA) using the Metagenomics 
RAST (MG-RAST) metagenomic analysis tool. R/c/ceffs/o-specific PCR was used for the characterization of rickettsial species. 
Almost full length of 16S rDNA was amplified in order to characterize unclassified bacterial sequences obtained in /. 
persulcatus female samples. The numbers of bacterial genera identified for the tick species were 71 (/. ovatus), 127 (/. 
persulcatus) and 59 (H. flava). Eighteen bacterial genera were commonly detected in all tick species. The predominant 
bacterial genus observed in all tick species was Coxiella. Spiroplasma was detected in Ixodes, and not in H. flava. PCA 
revealed that microbial populations in tick salivary glands were different between tick species, indicating that host 
specificities may play an important role in determining the microbial complement. Four female /. persulcatus samples 
contained a high abundance of several sequences belonging to Alphaproteobacteria symbionts. This study revealed the 
microbial populations within the salivary glands of three species of ticks, and the results will contribute to the knowledge 
and prediction of emerging tick-borne diseases. 
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Introduction 

Ticks (Acari: Ixodida) are globally one of the most important 
arthropod vectors of infectious diseases in animals and humans. 
They carry and transmit a range of bacterial, protozoan and viral 
pathogens, which are often zoonotic [1,2,3]. Opportunities for 
ticks to come into contact with humans and animals are increasing 
as their habitats change and their distribution widens. The 
incidence of tick-borne cases, including those emerging diseases, is 
on the rise [3,4]. Recent studies have identified microorganisms, 
such as Leptospira and Chlamydia, which had previously not been 
associated with ticks [5,6]. 

Ixodid ticks have four stages to their life cycle (egg, larva, 
nymph, and adult), and each post-embryotic phase requires blood 
to grow and molt. During feeding, ticks inject saliva into host 



animals to facilitate the uptake of blood. Tick saliva contains 
bioactive components that inhibit blood coagulation and the host 
immune system, and it also assists individuals to make a strong 
attachment to the skin of the host [7,8]. Tick-borne pathogens, 
such as Rickettsia, Ehrlichia and Anaplasma, concentrate within 
the salivary glands, and are transferred into the host animal during 
feeding on blood [9,10,11,12]. Electron microscopy has revealed 
the presence of non-pathogenic symbionts, such as Coxiella-like 
bacteria, in the salivary glands of ticks [13]. Symbionts have also 
been identified by molecular methods; however, the functions of 
these bacteria within tick salivary glands remain unclear 
[14,15,16,17,18]. Thus, the investigation of the microbial popu- 
lations in the tick salivary glands may lead to a better 
understanding of the microbes that could be transmitted to the 
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mammalian hosts together with tick saliva or that could play roles 
in tick biology such as by facilitating blood feeding. 

High throughput sequencing has provided an insight to the 
diversity of microbes associated with ticks [5,19,20,21]. The 
analysis of 16S ribosomal DNA (16S rDNA) amplicons by 
pyrosequencing is a high throughput technique that can be used 
to detect non-culturable microbes, and reveal entire populations of 
tick-borne microbes. The analysis of microbial populations found 
within tick salivary glands can provide important information to 
help predict emerging pathogens, assess potential risks, and 
understand the interactions between tick symbionts and patho- 
gens. 

In this study, tick samples were collected in the Shizuoka 
Prefecture, Japan, where Japanese spotted fever is endemic [22]. 
The aim of this work was to elucidate the microbial populations 
(including pathogens) associated with salivary glands of ticks using 
16S rDNA amplicon pyrosequencing technology. 

Materials and Methods 

Sample collection and DNA preparation 

Adult host-questing ticks of the species Ixodes ovatus, I. 
persulcatus, and Haemaphysalis flava were collected by flagging 
in mountainous areas of Shizuoka Prefecture from 2008 to 2010. 
No specific permissions were required for these locations and 
activities. Our field activities did not involve endangered or 
protected species. Table SI indicates information about the 
sampling sites. The sample numbers of /. ovatus, I. persulcatus, 
and H. flava used for this study were 24 (14 female, 10 male), 12 
(6 female, 6 male), and 5 (female only), respectively. Each tick was 
split into two parts (anterior and posterior) at the area between 
coxa 1 and coxa 2 using a microtome. The anterior part was then 
removed from the posterior part using sterile forceps. Since the 
salivary glands were attached with the anterior part, they could be 
removed from the tick carcass without damaging the midgut. The 
salivary glands were then collected into a sterile 1 .5 ml tube using 
sterile forceps, followed by washing with sterile PBS in order to 
minimize bacterial contamination. All dissection steps were 
performed under a stereomicroscope with great care to avoid 
the contamination of the midgut fluid. Genomic DNA was 
individually extracted using QJAamp DNA Mini kit (QIAGEN, 
Tokyo, Japan) according to the manufacturer's instructions, and 
stored at — 20°C. Samples of /. ovatus female and male, /. 
persulcatus female and male, and H. flava female are indicated by 
IOf, IOm, IPf, IPm, and HFf, respectively, throughout (Table 1). 

PCR amplification of V1 to V3 regions for 16S rDNA 
amplicon libraries 

The VI, V2, and V3 hyper-variable regions of bacterial 16S 
rDNA were amplified by PCR using the universal primers 27F (5'- 
X-AGAGTTTGATCMTGGCTCAG-3 ') and 518R (5'-AT- 
TACCGCGGCTGCTGG-3'), corresponding to positions 27 to 
518 of the Escherichia coli 16S rDNA [23,24]. The 27F primer 
contained ten bases of a multiplex identifier (MID) tag sequence 
denoted as 'X'. Primers 27F and 518R were modified with 5'- 
adapter A and 5' -adapter B sequences, respectively, for pyrose- 
quencing (Roche, Basel, Switzerland). PCR was performed in a 
total volume of 50 containing PCR Buffer, 0.2 |ll of Platinum 
Taq DNA polymerase (Life technologies, Tokyo, Japan), 0.2 mM 
of each primer, 1 (Xl of 10 mM dNTP mixture, 1.5 JLll of 50 mM 
MgCl 2 , and 1 |ll of template DNA. The PCR reaction was 
preceded with a 2 min denaturation step at 94°C, followed by 30 
cycles of 94°C for 30 sec, 55°C for 30 sec, and 72°C for 1 min. 
The quality of the PCR products (about 500 bp in length) were 



assessed by agarose (1%) gel electrophoresis, followed by 
purification using the Wizard SV Gel and PCR Clean-Up System 
(Promega, Tokyo, Japan). Concentration and quality of the 
amplicons were assessed with an Agilent 2100 BioAnalyzer 
(Agilent Technologies, Palo Alto, USA) using a DNA 1000 lab 
chip (Agilent Technologies). 

Pyrosequencing and data analysis 

Amplicon libraries with different MID tags were mixed and 
subjected to pyrosequencing using a 454 Genome Sequencer 
Junior (GS-Junior; Roche) following the manufacturer's protocol. 
The pyrosequencing data were deposited in DDBJ with accession 
no. DRA001731. The resulting data files (standard flowgram 
format, .sff files) were converted to FASTA files and sorted 
according to sample-specific MID tags using CLC Genomics 
Work Bench (CLC Bio, Tokyo, Japan). MID tag barcodes and 
primers were trimmed, then low-quality and short sequence reads 
(< 150 bp) were removed. DECIPHER'S Find Chimeras web tool 
(htqD://decipher.cee.wisc.edu/FindChimeras.html) was used to 
remove chimeric sequences. The remaining reads were phyloge- 
netically classified with the assistance of the Ribosomal Database 
Project (RDP) 16S Classifier version 10 (http://rdp.cme.msu.edu/ 
index.jsp), which can accurately and rapidly provide assignments 
for domains to the genus level. A comparative analysis of each 
sample was performed using the MG-RAST metagenomics 
analysis server employing the RDP dataset (lutp:// 
metagenomics.anl.gov/). Alpha diversity of each sample was also 
calculated using the MG-RAST server. Data sets were represented 
as the mean ± standard deviations (S.D.) after the Smirnov- 
Grubbs outlier test (a = 0.05). 

Conventional PCR methods 

Rickellsia-speciFic PCR amplification of the citrate synthase 
gene (gltA) using the primers RpCS877p (5'- 
GGGGGCCTGCTCACGGCGG-3') and RpCS1273r (5'-CAT- 
AAC C AGTGTAAAGCTG-3 ') [25] was performed on 22 sam- 
ples that were highlighted by RPD analysis as containing the genus 
Rickettsia. PCR was performed in a total volume of 25 |ll 
containing PCR Buffer for KOD-Plus-Neo, 0.5 |il of KOD-Plus- 
Neo DNA polymerase (Toyobo, Tokyo, Japan), 0.3 mM of each 
primer, 2.5 nj of 2 mM dNTP mixture, 1.5 |il of 25 mM MgS0 4 , 
and 1 (4.1 of template DNA. The reaction started with a 
denaturation step at 94°C for 2 min, followed by 40 cycles of 
94°C for 15 sec, 54°C for 30 sec, and 68°C for 30 sec, and a final 
extension step at 68°C for 2 min. PCR products were purified 
using ExoSap-IT (Affymetrix, Tokyo, Japan) according to the 
manufacturer's instructions. Cycle sequencing was performed 
using BigDye v3.1 terminator chemistry (Applied Biosystems, 
Tokyo, Japan) and the same primers. Sequencing products were 
assessed using a 3130x1 Genetic Analyzer (Life Technologies, 
Tokyo, Japan). 

RDP analysis did not classify all sequences to the genus level. To 
characterize these sequences, the near full-length 1 6S rDNA gene 
was obtained from four /. persulcatus female samples by PCR 
using the universal primers fD 1 (5'-AGAGTTTGATCCTGGCT- 
CAG-3') and Rp2 (5 ' - AC GGC T AC C TTGTT AC G AC TT- 3 ' ) 
[26]. PCR was carried out in a final volume of 50 |ll PCR Buffer 
for KOD-Plus-Neo, 1 \l\ of KOD-Plus-Neo DNA polymerase, 
0.3 mM of each primer, 5 pi of 2 mM dNTP mixture, 3 jal of 
25 mM MgS0 4 , and 2 |ll of template DNA. PCR conditions 
started with a denaturation step at 94°C for 2 min, followed by 40 
cycles of 94°C for 15 sec, 55°C for 30 sec, and 68°C for 45 sec, 
and a final extension step at 68°C for 2 min. Quality of the PCR 
products (approx. 1400 bp) was assessed by agarose (1%) gel 
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Table 1. 


Sequence results and number of detected genera. 








Sample ID 


Tick species 


Sex 


No. of sequence reads 


No. of genera 


lOfl 


/. ovatus 


female 


5,664 


10 


IOf2 


1. ovatus 


female 


4,484 


10 


IOf3 


1. ovatus 


female 


3,498 


6 


IOf4 


1. ovatus 


female 


4,712 


7 


IOf5 


1. ovatus 


female 


5,591 


3 


IOf6 


1. ovatus 


female 


4,200 


5 


IOf7 


1. ovatus 


female 


5,030 


6 


IOf8 


1. ovatus 


female 


5,634 


23 


IOf9 


1. ovatus 


female 


7,643 


8 


lOflO 


1. ovatus 


female 


5,636 


7 


lOfll 


1. ovatus 


female 


4,049 


2 


IOf12 


1. ovatus 


female 


7,275 


14 


IOf13 


1. ovatus 


female 


3,351 


3 


IOf14 


1. ovatus 


female 


7,048 


7 


10ml 


1. ovatus 


male 


4,986 


13 


IOm2 


1. ovatus 


male 


3,790 


12 


IOm3 


1. ovatus 


male 


7,916 


22 


IOm4 


1. ovatus 


male 


3,844 


5 


IOm5 


1. ovatus 


male 


6,340 


18 


IOm6 


1. ovatus 


male 


7,130 


22 


IOm7 


1. ovatus 


male 


6,176 


22 


IOm8 


1. ovatus 


male 


9,788 


28 


IOm9 


1. ovatus 


male 


9,628 


12 


lOmlO 


1. ovatus 


male 


7,170 


17 


IPfl 


1. persulcatus 


female 


8,964 


38 


IPf2 


1. persulcatus 


female 


3,599 


16 


IPf3 


1. persulcatus 


female 


7,085 


42 


IPf4 


1. persulcatus 


female 


8,242 


40 


IPf5 


1. persulcatus 


female 


7,943 


25 


IPf6 


1. persulcatus 


female 


1 0,506 


18 


IPml 


1. persulcatus 


male 


7,414 


55 


IPm2 


1. persulcatus 


male 


8,173 


19 


IPm3 


1. persulcatus 


male 


1 0,803 


26 


IPm4 


1. persulcatus 


male 


10,144 


34 


IPmS 


1. persulcatus 


male 


16,117 


29 


IPm6 


1. persulcatus 


male 


9,221 


40 


HFfl 


H. flava 


female 


6,438 


3 


HFf2 


H. flava 


female 


8,339 


44 


HFf3 


H. flava 


female 


6,204 


5 


HFf4 


H. flava 


female 


10,017 


18 


HFf5 


H. flava 


female 


8,294 


24 


doi:1 0.1 371 /journal.pone.01 03961 .t001 



electrophoresis, followed by purification of the products using the Sanger sequencing data analysis 

Wizard SV Gel and PCR Clean-Up System (Promega). PCR Sanger sequencing data were analyzed using GENETYX 

products were A-tailed and then cloned with TA cloning plasmids version 9.1 (GENETYX Corporation, Tokyo, Japan). The 
pGEM-T Easy (Promega). Ten clones per sample were randomly GenBank accession numbers for the gltA sequences are 
selected and sequenced. AB911107 to AB911109, and the 16S rDNA sequences 

AB906824 to AB906829. Sequences were compared with those 
in public databases using nucleotide BLAST at NCBI website 
(http://blast.ncbi.nlm.nih.gov/Blast.cgi). 
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Figure 1. Relative abundances of different bacterial genera in the salivary glands of (A) /. ovatus, (B) /. persulcatus and (C) H. flava. All 
genera with less than 1.0% contribution were pooled into one group and labelled "others". 
doi:1 0.1 371 /journal.pone.01 03961 .g001 
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Figure 2. Venn diagram of all 163 identified genera distributed across the tick species and sex. 

doi:1 0.1 371 /journal.pone.01 03961 .g002 



Phylogenetic analysis was conducted using MEGA version 6.05 
[27]. The universal 16S rDNA sequences were aligned with those 
of closely related bacteria in GenBank using ClustalW and a 
maximum likelihood phylogram was constructed. 

Results 

Classification and quantification of bacterial taxa 

Between 3,351 and 9,788 sequences were obtained for 
individual /. ovatus, of which almost 98% were assigned to the 
genus level (Table 1) (Figure 1A). A total of 71 bacterial genera 
were detected in /. ovatus, with 59 found in males and 37 in 
females. The two dominant bacterial genera were Spiroplasma and 
Coxiella, and these accounted for more than 90% of the bacterial 
community in ticks, except for a single /. ovatus female and 3 /. 
ovatus males (Figure 1A). Rickettsia (genus contains known tick- 
borne pathogens R. japonica and R. helvetica) was detected in ten 
samples, and Ehrlichia (genus contains known tick-borne patho- 
gens E. chaffeensis and E. muris) was detected in two samples. 

Between 3,599 and 16,117 sequences were recorded for 
individual /. persulcatus, with almost 82% assigned to the genus 
level (Table 1) (Figure IB), except for those of four /. persulcatus 
females. Over 80% of the reads in these samples were unclassified 
at the genus level. At the phylum level, these reads were classified 
as Proteobacteria or Alphaproteobacteria by the RDP classifier 
(data not shown). There were 127 different bacterial genera 
detected in /. persulcatus, of which 92 were detected in males, and 



81 in females. Rickettsia was detected in nine /. persulcatus (4 
female, 5 male) individuals and Ehrlichia was detected in a single 
/. persulcatus male (IPm5). 

Between 6,204 and 10,017 sequences were obtained for 
individual H. flava, of which almost 97% were identified to the 
genus level (Table 1) (Figure 1C). A total of 59 different bacterial 
genera were detected, and Coxiella accounted for more than 90% 
of the microbial population in all samples. Spiroplasma was not 
detected in any individuals of H. flava, despite appearing in all 
Ixodes samples. Rickettsia spp. were detected in three H. flava 
females, and no sample contained Ehrlichia spp. The details of 
microbial population analysis for each sample are presented in 
Table S2. 

A summarized diagram of the number of bacterial genera 
detected in each tick group is presented in Figure 2. Out of 163 
different genera identified, 18 were detected in all tick groups. 
These were Acinetobacterium, Arcicella, Burkholdelia, Corynebac- 
terium, Coxiella, Cryobacter, Curvibacterium, Flavobacterium, 
Limnohabitas, Methylobacterium, Novosphingomonas, Polynucleo- 
bacter, Propionilbacterium, Pseudomonas, Rickettsia, Sphingomo- 
nas, Staphylococcus, and Sterptophyta. Some bacterial genera were 
uniquely associated with tick species or sex, i.e., IOf (1 genus), 
IOm (13 genera), IPf (24), IPm (35) and HFf (19). 
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Figure 3. Principal component analysis of the bacterial composition in each tick sample. The plots were generated using the MG-RAST 
server. Each tick sample is shown in a different color depending on the species and sex of the tick; lOf, lOm, IPf, IPm, and HFf are respectively, shown 
in red, green, blue, purple, and yellow. The plots derived from the same tick species are highlighted in circles; /. ovatus (IO), /. persulcatus (IP), and H. 
flava (HF) are, respectively, highlighted in red, blue, and yellow circles. 
doi:1 0.1 371 /journal.pone.01 03961 .g003 



Comparison of microbiomes in salivary glands between 
tick species 

PCA was performed using the MG-RAST server with 
normalized values and Bray-Curtis distance (Figure 3) for each 
tick sample. The microbial community composition of each 
sample clustered approximately according to tick species. The 
microbial populations of Ixodes and Haemaphysalis were com- 
pletely separated by PC02. The microbial community composi- 
tion of Haemaphysalis ticks was broadly distributed along PCOl; 
however, in /. ovatus and /. persulcatus microbial populations 
were more distinct, but with some overlap within this component. 

Alpha diversity for each sample was calculated using the MG- 
RAST server (Figure 4). Smirnov-Grubbs's outlier test (ot = 0.05) 
was used before the calculation of means and S.D. IPml was 
identified as an oudier and removed in the calculation for the 
mean value of IPm alpha diversity. Mean values were 5.75 ±1.19 
(IOf), 5.33±0.72 (IOm), 4.97±1.25 (IPf), 3.11±0.55 (IPm) and 
2.14±0.32 (HFf). 

Sequencing of gltA 

The gltA gene was detected in 1 1 out of 22 samples previously 
identified as containing the genus Rickettsia. Samples that were 
gltA -positive tended to have a greater abundance of rickettsial 
bacteria than those that were negative (Figure 5). All gltA -positive 
samples were subjected to sequencing analysis. Each gltA -positive 



sample contained only one sequence type, indicating that 
individual ticks harbored bacteria carrying a single gltA allele. 
From 1 1 tick samples, three different gltA sequences were 
identified, and BLAST searches showed the highest identities 
(99.8% to 100%) with R. asiatica, R. helvetica, and uncultured 
Rickettsia sp. (Table 2). 

Sequencing of unclassified bacterial 16S rDNA 

PCR fragments (1400 bp) were generated using universal 
primers to resolve the identities of sequences detected in four 
individuals of /. persulcatus. Between six and nine clones per 
sample were classified into Alphaproteobacteria (data not shown) 
based on BLASTn similarity searches. All the clones analyzed 
from two individual ticks were the same sequence type. There 
were four different sequence types in one individual, and two in 
another. These showed highest identities (99.5% to 99.7%) with 
uncultured Rickettsiales previously reported from /. persulcatus 
(GenBank accession number AF497583). 

Molecular phylogenetic analysis revealed that the Alphaproteo- 
bacteria from four /. persulcatus females clustered together within 
a single clade. This clade contains Candidatus Lariskella 
arthropodarum identified in several stinkbug species [Arocatus 
melanostomus, Nysius plebeius, and Physopelta gutta) and Rickett- 
siales derived from flea (Xenopsylla cheopis) and ticks (/. ovatus and 
/. persulcatus) [28,29,30,31] (Figure 6). 
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Figure 4. Alpha diversity calculated for each tick sample. The alpha diversity of each tick sample was calculated using the MG-RAST server. 
The mean value obtained for each tick group is represented by the horizontal line. Mean alpha diversity values: IOf (5.75), lOm (5.33), IPf (4.97), IPm 
(3.11), and HFf (2.14). 
doi:1 0.1 371 /journal.pone.01 03961 .g004 



Discussion 

The aim of this study was to assess and compare the diversity of 
bacterial populations within the salivary glands of /. ovatus, I. 
persulcatus, and H. flava. This is the first reported study of 
bacterial populations found in the salivary glands of different tick 
species. This metagenomic approach revealed bacterial popula- 
tions totaling 163 different genera found in tick salivary glands. 
These included the genera of tick-borne pathogens such as 
Ehrlichia and Rickettsia. Further identification using species- 
specific PCR would be needed to clarify the presence of the tick- 
borne bacterial pathogens, such as E. muris, E. chaffeensis, R. 
japonica and R. helvetica, in the ticks used in this study 
[32,33,34,35]. This combination of detection approaches may be 
useful for the screening and detection of possible pathogens in 
arthropod vectors. 

Rickettsia was detected in 22 of the 41 (53.6%) samples by 16S 
rDNA amplicon pyrosequencing; however, only half of the 22 
positive samples were positive with gllA PCR. This may be 
attributed to the relative amounts of rickettsial DNA in the PCR 
templates, where gltA PCR-positive samples tended to contain a 
higher proportion of rickettsial DNA than those that were negative 
(Figure 5). However, there were two gltA -negative samples (IPml 
and IPf4) that had higher proportions of rickettsial DNA than a 
gllA -positive sample (IOf4). This suggests that the sensitivity of 



conventional gllA PCR may be affected by the other factors such 
as the resolving power of agarose gel electrophoresis and the 
presence of PCR inhibitory components in samples [36,37]. We 
suggest that a 16S rDNA amplicon pyrosequencing approach is a 
more sensitive method to detect specific pathogens. 

Analysis of the gltA gene sequences from /. ovatus and /. 
persulcatus revealed that they belonged to R. asiatica and R. 
helvetica, respectively (Table 2). This result agrees with previous 
findings about the potential of the ticks to act as vectors for these 
rickettsia in Japan [38]. R. helvetica belongs to the spotted fever 
group of rickettsia and is a causative agent of febrile illness. A 
human case associated with this pathogens has been reported 
elsewhere [39,40]. There was a high abundance (>70%) of this 
rickettsial species in some /. persulcatus samples (Figure IB), 
suggesting that it is well adapted to the salivary glands of ticks, and 
waiting for transmission to mammalian hosts. In addition to 
pathogenic strains, the genus Rickettsia also contains symbionts 
associated with ticks. Rickettsia-like symbionts can influence the 
tick physiology, population dynamics, and the transmission of 
other pathogenic Rickettsia spp [4 1 ,42] . 

Coxiella burnetii and Coxiella-Vike endosymbionts have been 
identified in several tick genera, including Dermacentor, Ixodes, 
Haemaphysalis and Rhipicephalus [43,44,45,46]. Coxiella-like 
endosymbionts have been located at high densities in the salivary 
glands of the lone star tick (Amblyomma americanum) using 
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Figure 5. Comparison of the relative abundance of rickettsial sequences estimated by 1 6S amplicon analysis and the results of gltA 
PCR. Vertical axis represents the relative abundance of rickettsial sequences calculated from the data obtained from 1 6S amplicon analysis. Blue dots 
represent samples in which Rickettsia was detected by both 16S amplicon analysis and gltA PCR. Red dots represent samples in which Rickettsia was 
detected by 16S amplicon analysis but not by gltA PCR. The plots with relative abundance values between 0% and 5% are shown in the magnified 
graph provided in the right column. 
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fluorescence in situ hybridization [13]. The findings in this study three species of tick. The dominant presence of Coxiella in the 
also highlighted the presence of Coxiella in the salivary glands of salivary glands of ticks warrants further investigation to resolve 



Table 2. Summary of gltA sequencing. 



Identity with 
reference (no. 



Sequence 


Tick sample 


matched/no. 


Reference 


Rickettsia 


GenBank 


ID 


ID 


nucleotides) 


GenBank no. 


species 


no. 


gltAJOfl 


lOfl 


99.8% (438/439) 


AB297808 


R. asiatica 


AB911107 


gltA_IOf3 


IOf3 


99.8% (438/439) 


AB297808 


R. asiatica 


AB911107 


gltA_IOf4 


IOf4 


99.8% (438/439) 


AB297808 


R. asiatica 


AB911107 


gltA_IOm2 


IOm2 


99.8% (438/439) 


AF394901 


R. asiatica 


AB911107 


gltA_IOm3 


IOm3 


99.8% (438/439) 


AF394901 


R. asiatica 


AB911107 


gltA_IOm4 


IOm4 


99.8% (438/439) 


AF394901 


R. asiatica 


AB911107 


gltAJPfl 


IPfl 


99.8% (438/439) 


U59723 


R. helvetica 


AB911108 


gltA_IPf6 


IPf6 


1 00% (394/394) 


JN849396 


Uncultured 


AB911109 










Rickettsia sp. 




gltA_IPm4 


IPm4 


99.8% (438/439) 


U59723 


R. helvetica 


AB911108 


gltA_IPm5 


IPm5 


99.8% (438/439) 


U59723 


R. helvetica 


AB911108 


gltA_IPm6 


IPm6 


99.8% (438/439) 


U59723 


R. helvetica 


AB911108 
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Figure 6. Phylogenetic analysis of the 16S rDNA sequences of unclassified bacteria from IPf2, IPf3, IPf4, and IPf5 using maximum 
likelihood method. The tree is rooted with the Escherichia coli. All bootstrap values from 1000 replications are shown on interior branch nodes. 
doi:1 0.1 371 /journal.pone.01 03961 .g006 



their potential roles in tick biology, particularly blood-sucking 
behavior, and their interaction with other microbes. 

The genus Spiroplasma contains a wide diversity of often 
unnamed or poorly characterized species, including non-patho- 
genic, symbiotic, and pathogenic organisms associated with a wide 
variety of arthropods. Symbiotic Spiroplasma has a close 
association with, and can affect the behavior of, their host 
arthropods. For example, Hurst et al. (2000) reported the 
preferential killing of males by Spiroplasma; when female insects 
(e.g., the butterfly Danaus chrysippus) are infected, the broods are 
female-biased because the infected male progeny die during 
embryogenesis [47]. One Spiroplasma sp. has been reported in 
ticks [48,49], and it has also been associated with transmissible 
spongiform encephalopathy in humans and ruminants, although 
its role in the pathology of the host has not been clarified [50] . In 
this study, Spiroplasma was detected in Ixodes ticks, and not in H. 
flava (Figure 1). Previous research reported the genera Spiro- 
plasma and the closely related Mycoplasma in several tick species 
in Japan [51]. The pathogenicity of Spiroplasma harbored in ticks 
in Japan is not known yet. 

Results from the PCA of sequences indicated that microbial 
population structures in the salivary glands of ticks were different, 
and that samples from the same species of tick clustered together 
(Figure 3). Ticks can acquire microorganisms through a variety of 
ways, such as transovarial transmission, and from the environ- 
ment, host animals during blood feeding, and mating partners. For 
microorganisms to exist in the salivary glands, they need to 



migrate from the midgut and enter the glands. The establishment 
of microorganisms within ticks can depend on the interactions 
between particular microbes, ticks and other symbioses [41,52,53]. 
The differences in the microbial populations within the salivary 
glands of tick species in this study were attributed to these 
complicated factors. 

Previous studies revealed that tick microbial populations were 
different between developmental stages (egg, nymph, and adult) 
[19,21]. The bacterial compositions also differed between organs, 
such as between midgut and ovary [19]. Some bacterial species, 
for instance Borrelia burgdorferi that is a causative agent of Lyme 
disease, exist in the midgut of the tick, moving into the salivary 
glands when stimulated by feeding on blood [54,55]. For better 
understanding of microbial interactions with ticks as well as the 
potential pathogens transmitted by ticks, further study should 
include the comparison of the microbes between salivary glands 
and other organs. The analysis of the dynamics of microbial 
community composition during the process of feeding on blood 
may also uncover the roles of tick microbes. 

The mean alpha diversity value (Figure 4) was greater for the 
female /. ovatus (5.61) than that of male (5.31). This rank order 
was also recorded for female (5.02) and male (3.38) /. persulcatus 
ticks, and may imply that some bacterial species preferentially 
select the gender of ticks. There may be some strategic biological 
relevance in the transmission of bacteria to mammalian hosts 
because female ticks feed for a longer period of time than males. 
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The total number of bacterial genera (Table 1) detected in /. 
persulcatus (127) was greater than in /. ovatus (71). 

Several /. persulcatus females contained unclassified bacteria 
belonging to the Proteobacteria and Alphaproteobacteria (Fig- 
ure IB). Based on the analysis of the nearly complete 16S rDNA 
sequences, the unclassified bacterial were classified into a single 
phylogenetic clade, which was recently proposed as a "Candidatus 
L. arthropodarum" clade [30]. This clade also includes Rickett- 
siales bacterium previously found in blood and biopsy samples of 
the patients with an acute fever disease, etiologically linked with 
tick bites [31]. The relationships between these microorganisms 
and their arthropod hosts are not clear, and their potential to act 
as causative agents of emerging tick-borne mammalian diseases 
warrants further investigation. 
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